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Superconductive  quantum  interference-devices  (SQUIDs)  offer  new  technology  for 
locating  materials  flaws  electromagnetically  that  promises  to  increase  sensitivity, 
depth  of  magnetic  flux  enables  use  of  microscopic  pickup  loops  in  a  gradiometer 
configuration  to  give  high  resolution.  A  cryogenic  umbilical  connects  pickup  loops 
to  a  remote  cryostat  housing  SQUID  sensors  to  ease  scanning.  A  pair  of  drive  coils 
a  few  millimeters  in  radius  that  encircle  pickup  loops  forming  a  coplanar 
gradiometer  1  mm  or  less  in  radius  comprise  a  superconductive  microprobe.  It 
provides  a  depth  of  field  of  several  millimeters  to  a  0.1  mm  flaw  in  an  aluminum 
plate,  when  operating  with  a  drive  current  a  1  A  oscillating  at  a  frequency  of 
1kHz.  Its  field  of  view  ranges  to  several  millimeters,  for  flaws  a  few  millimeters 
deep,  and  its  horizontal  resolution  is  1  mm  or  so,  for  flaw  depths  out  to  its  depth 
of  field.  An  array  of  microprobes  form  receptors  much  like  rods  in  the  retina  of  a 
|  magnetic  eye.  The  eye  leads  to  an  electromagnetic  microscope  for  imaging  internal 
flaws  in  aluminum  plates.  It  gives  multiple  images  that  enable  resolving  depth  of 
a  0.1  mm  flaw  to  a  few  tenths  of  a  millimeter  with  a  horizontal  resolution  of  one 
millimeter  or  so.  f i< p  ) 
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SUMMARY 


Superconductive  quantum  interference  devices  (SQUIDs)  offer  new  technology  for 
locating  material  flaws  electromagnetically  that  promises  to  increase  sensitivity,  depth  of 
field,  and  field  of  view  and  to  enhance  resolution.  The  ultrahigh  sensitivity  of  SQUIDs  to 
magnetic  flux  enables  use  of  microscopic  pickup  loops  in  a  gradiometer  configuration  to 
give  high  resolution.  A  cryogenic  umbilical  connects  pickup  loops  to  a  remote  cryostat 
nousing  SQUID  sensors  to  ease  scanning. 

A  pair  of  drive  coils  a  few  millimeters  in  radius  that  encircle  pickup  loops  forming  a 
coplanar  gradiometer  1  mm  or  less  in  radius  comprise  a  superconductive  microprobe.  It 
provides  a  depth  of  field  of  several  millimeters  to  a  0.1  mm  flaw  in  an  aluminum  plate, 
when  operating  with  a  drive  current  of  1  A  oscillating  at  a  frequency  of  1  kHz.  Its  field 
of  view  ranges  to  several  millimeters,  for  Haws  a  few  millimeters  deep,  and  its  horizontal 
resolution  is  1  mm  or  so,  for  flaw  depths  out  to  its  depth  of  field. 

An  array  of  microprobes  form  receptors  much  like  rods  in  the  retina  of  a  magnetic 
eye.  The  eye  leads  to  an  electromagnetic  microscope  for  imaging  internal  flaws  in 
aluminum  plates.  It  gives  multiple  images  that  enable  resolving  depth  of  a  0.1  mm  flaw  to 
a  few  tenths  of  a  millimeter  with  a  horizontal  resolution  of  one  millimeter  or  so. 
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1.0  INTRODUCTION 


Superconductive  quantum  interference  devices  (SQUIDs)  offer  new  technology  for 
locating  material  flaws  electromagnetically  that  promises  to  increase  sensitivity  and  depth 
of  field,  to  enhance  resolution,  and  to  enable  multivariate  measurements  with  a  single 
system.  They  can  measure  changes  in  magnetic  flux  as  small  as  10 '4  flux  quanta  (a  flux 
quantum  is  2.07  X  10  "15  Tm*)  at  frequencies  from  0.1  Hz  to  several  kilohertz.  A 
conventional  pickup  coil,  operating  at  5  MHz  with  a  sensitivity  of  1  fiV,  requires  about  105 
turns  of  wire  to  give  an  equal  flux  sensitivity,  or  alternatively  a  100  turn  coil  must  operate 
at  a  frequency  of  5  GHz,  which  is  equally  impractical.  A  SQUID,  however,  needs  but  a 
single  loop  one  millimeter  in  diameter. 

The  ultrahigh  sensitivity  of  SQUIDs  to  magnetic  flux  allows  use  of  microscopic 
pickup  loops  in  a  gradiometer  configuration  to  give  high  resolution,  arrays  of  closely 
spaced  loops  can  provide  high  scan  rates  and  image  flawed  regions.  High  sensitivity  also 
allows  operating  at  standoff  distances  of  several  millimeters,  which  increases  aperture  and 
hence  scan  rate.  Surveys  can  be  made  rapidly  at  large  standoff  and  low  resolution  followed 
by  mapping  flawed  regions  at  small  standoff  and  high  resolution. 

Because  their  sensitivity  is  independent  of  frequency  over  a  wide  range,  SQUIDs 
enable  using  low  frequencies  to  increase  depth  of  field  as  well  as  operating  at  multiple 
frequencies  to  scan  over  depth.  They  also  can  measure  minute  magnetic  inhomogeneities  of 
a  material  and  so  evaluate  magnetic  properties  such  as  remanent  magnetization,  magnetic 
permeability,  coercive  field  strength,  and  Barkhausen  noise.  The  full  range  of 
measurements  of  magnetic  and  electrical  properties  available  provides  the  versatility  of 
multivariate  measurements  with  a  single  system. 

Phase  I  research  reported  here  uses  model  analyses  to  quantify  advances  anticipated 
from  using  SQUID  technology  for  nondestructive  evaluation  of  materials.  It  demonstrates 
feasiblity  of  the  design  of  an  eddy  current  probe  that  uses  microscopic,  coplanar  pickup 
loops  in  a  gradiometer  configuraton  to  attain  a  resolution  of  one  millimeter  and  a  depth  of 
field  of  several  millimeters  to  a  0.1  mm  flaw  in  an  aluminum  plate,  at  operating  frequencies 
from  a  few  hundred  to  a  few  thousand  Hertz.  The  design  leads  to  an  array  of  closely 
spaced  microprobes  suited  to  scanning  plates  and  tubes  rapidly  and  with  high  resolution. 

To  realize  the  concept  of  a  microscopic  probe,  we  set  forth  an  innovative  cryogenic 
design  that  attaches  the  probe  to  a  compact  cryostat,  housing  SQUID  sensors,  through  an 
umbilical  formed  by  a  cryogenic  heat  pipe.  A  cluster  of  microprobes  forms  an  array  much 
like  rods  in  the  retina  of  a  magnetic  eye.  The  eye  leads  to  an  electromagnetic  microscope 
the  gives  multiple  images  of  internal  flaws,  allowing  resolution  in  depth  of  a  few  tenths  of  a 
millimeter  to  a  0.1  mm  flaw  and  a  horizontal  resolution  of  1  mm  or  so. 
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2.0  CONFIGURATIONS  OF  A  MICROPROBE 


Figure  2.1  shows  a  configuration  of  source  coils  and  pickup  loops  that  comprise  a 
microprobe.  They  form  a  coplanar  nest  of  concentric  circles.  The  two  outermost  coils 
carry  an  oscillating  electric  current  to  drive  eddy  currents.  Current  in  the  inner  drive  coil 
opposes  current  in  the  outer  coil  in  order  to  suppress  interference  from  the  source  at  the 
receiver,  formed  by  the  inner  most  loops.  Loops  of  the  receiver  wind  oppositely  to  form  a 
planar  gradiometer  of  order  n,  where  n  is  the  number  of  loops.  Figure  2.2  shows  receiver 
windings  divided  into  pairs  of  opposing  quadrants  to  form  a  biaxial  gradiometer  of  order  n. 

We  first  specify  geometry  of  receiver  loops  forming  a  gradiometer  of  order  n  and 
then  give  specifications  of  source  coils. 

2.1  RECEIVER  LOOPS 


Receiver  loops  form  a  planar  gradiometer  of  order  n,  where  n  is  the  number  of 
pickup  loops.  A  loop  has  Ni  turns  and  radius  For  axial  symmetry,  the  vector 

potential,  A,  of  a  magnetic  flux  density,  B  =  ?  x  A  ,  threading  the  loops  points  in  the 

direction  of  cylindrical  coordinates  r,  <J>,  and  z,  where  z  points  along  the  common  axis 
of  the  loops.  The  expression 

A*,(t)  =  LCn  r"  (21a) 


expands  magnitude  of  the  vector  potential  in  powers  of  the  radial  coordinate  r.  For 


example,  a  magnetic  dipole,  Mz,  at  position  z  along  the  z  axis  gives  an  expansion  expressed 
by  the  relation 

00 


(2.1b) 


n  »  0 


where  n0  =  4tt  *  10*7  H/m  and  double  factorials  denote  products  of  odd  and  even  integers. 

The  magnetic  flux  threading  an  ith  loop  is  2tti  Ni  A^(rj),  so  Equation  2.1a  says 
that  the  relation 


R»=  Y  Ni  2ir  S  Cn  (ri)n+1 

mmt  n  =  l 

i  *  i 


(2.2a) 


gives  the  net  flux,  RB  ,  threading  a  receiver  formed  by  m  concentric  loops  with  n_i  >  ri 
>  0.  The  net  flux  or  response  Rffl  vanishes  to  order  rm  for  a  gradiometer  of  m th  order, 
so  Equation  2.2a  tells  us  that  the  relations 


^  Ni(ri)n+1  =  0  ,  for  n  =  1,  2,  3,..,  m-1,  (2.2b) 

i  *  1 

specify  ratios  of  radii  and  number  of  turns  of  its  m  concentric  loops. 
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Figure  2.1 

A  coplanar  nest  of  drive  coils  and  pickup  loops  forming  a  planar  gradiometer  of 
order  n.  The  two  outermost  loops,  numbered  0  and  1,  are  drive  coils  that  carry  an 
oscillating  current  to  excite  eddy  currents.  Their  currents  flow  oppositely,  as  marked  by 
arrows,  in  order  to  null  interference  from  drive  currents  at  the  receiver,  formed  by  n 
innermost  loops,  numbered  2  through  2+n.  Arrows  show  that  they  wind  oppositely  in 
order  to  give  a  null  response  to  magnetic  fields  radially  uniform  to  order  n-2.  Triple  dots, 
•  •  •  ,  between  receiver  loops  denote  continuation  of  windings  to  loop  n-i-2  . 


Figure  2.2 

Receiver  pickup  loops  forming  a  planar,  biaxial  gradiometer  of  order  n.  Opposite 
quadrants  form  pickup  loops  of  an  axis,  as  marked  by  heavy  lines.  Dashed  lines  mark 
pickup  loops  of  the  second  axis.  Arrows  show  winding  of  adjacent  loops  is  opposite. 
Closures  for  each  loop  run  parallel  along  diameters,  and  crossovers  are  at  the  center. 
Triple  dots,  •  •  •  ,  between  loops  denote  continuation  of  windings  to  loop  n+2  . 
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2.1.2  Second  Order  G radiometer 

For  a  gradiometer  of  second  order,  m  =  2,  and  Equation  2.2b  gives  the  relation 


Nir2  +  N2  r2  =  0 


(2.3) 


that  specifies  the  ratio  of  radii  and  number  of  turns  of  the  two  loops.  Namely, 

N2  r3  =  —1,  for  ri  =  1  and  Ni  =  1.  (Number  of  turns  is  positive  for  a  right  handed  winding 

and  negative  for  a  left  handed  winding.)  Equation  2.1b  then  tells  us  that  response  of  a 
second  order  gradiometer  to  a  dipole  at  position  z  on  the  axis  is  proportional  to  z'5  ,  to 
lowest  order. 

2.1.3  Third  Order  Gradiometer 

For  a  gradiometer  of  third  order,  m  =  3,  and  the  relations 

1  +  N2  i\  +  N3  rf  =  0  (2.4a) 

and 

1  +  N2  r*  +  N3  t\  =  0  (2.4b) 

constrain  number  of  turns  and  radii  of  the  two  innermost  loops,  with  Ni  =  1  and  rt  =  1. 
For  N2  =  —2  and  N3  =  2,  Equations  2.4a  and  2.4b  require  r2  =  0.8405  and  r3  =  0.4542. 

Response  to  a  dipole  on  the  axis  is  again  proportional  to  z*5  ,  to  lowest  order,  for  a 
gradiometer  of  third  order,  because  the  response  to  a  dipole  contains  even  powers  of  r 
alone. 


2.1.4  Fourth  Order  Gradiometer 

For  a  gradiometer  of  fourth  order,  m  =  4,  and  Equation  2.2b  gives  the  relations 

1  +  N2  rj  +  N3  r3  +  N4  r|  =  0  ,  (2.5a) 

1  +  N2  rj  +  N3  r3  +  N4  rj  =  0  ,  (2.5b) 

and 

1  +  N2  i\  +  N3  rj  +  N4  r{  =  0  ,  (2.5c) 

that  constrain  number  of  turns  and  radii  of  the  three  innermost  loops,  for  Nj  =  1  and 
r  1  =  1.  For  N2  =  -2,  N3  =  2,  and  N4  =  -2,  they  require  r2  =  0.899253,  r3  =  0.635352,  and 
r4  =  0.308249. 

Equation  2.1b  tells  us  that  response  to  a  dipole  on  the  axis  is  proportional  to  z-7,  to 
lowest  order,  for  a  gradiometer  of  fourth  order. 
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2.2  SOURCE  COILS 


The  source  comprises  two  coplanar  coils  that  encircle  the  receiver.  Electric  current 
oscillating  in  the  inner  coil  opposes  current  in  the  outer  coil  in  order  to  suppress 
interference  from  the  source  at  the  receiver.  Current,  Ii ,  in  the  inner  coil  balances  current 
in  the  outer  coil,  I0  ,  so  interference  from  the  source  vanishes.  We  set  the  current  so  that 
Ii  =  -0lo  ,  where  0  is  a  balance  factor. 

The  balance  factor  0  depends  on  mutual  inductances,  Mij(h)  ,  between  source 
loops,  numbered  0  and  1,  and  receiver  loops,  numbered  2  through  m  +  2,  as  expressed  by 
the  relation 


a  _  Mo2(h)  +  Mo3(h)  +  *  *  *  +  Mon 
P~  M12(h)  +  Mij(h)  +  •••  +  Mm  ’ 

with  n  =  m  +  2.  Here,  the  expression  (Appendix  B,  Equation  B4.a) 


My(h)  =  (loir  Ntfi  Njrj 


/  [ 1  +  [l?T^je'"2kh  ]  J'<krt>  dk 


(2.6a) 


(2.6b) 


gives  the  mutual  inductance  between  coplanar,  concentric  loops  of  radii  ri  and  rj  at 
height  h  above  an  electrically  conducting  plate,  where  ?2  =  k2  +  i (i0<rv  ,  a  is  electrical 
conductivity  of  the  plate,  u  is  the  carrier  frequency,  po  =  4x  *  10‘7  H/m  ,  and  Ji(x)  is  a 
Bessel  function  of  the  first  kind  and  order  one. 

We  choose  the  number  of  turns  and  radius,  Nt  and  ri  ,  of  the  inner  drive  coil  so 
that  0  is  unity  far  above  the  plate  (h  >  0)  and/or  at  low  frequency  (7  »  k),  where  the 
exponential  term  in  Equation  2.6b  vanishes.  The  expression  (Appendix  B,  Equation  B4.c) 

lim  Mij  =  M99  =  no  Ni  Nj  [  (rf+rj)  K(p9.  )  -  (ri+rjp  E(p?.  )  1  (2.7a) 

then  gives  mutual  inductances,  where  K(p)  and  E(p)  are  complete  elliptic  integrals  of  the 

first  and  second  kinds,  respectively,  and  p°i}  =  2  /rirj/(ri  +  rj).  Near  the  plate,  amplitude 

and  phase  of  the  factor  0  specify  current  in  the  inner  drive  coil  needed  to  null  interference 
at  the  receiver. 
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3.0  FLAW  SIGNAL  DESCRIPTION 


In  order  to  examine  response  of  a  microprobe  to  a  flaw,  we  first  determine  response 
of  a  single  receiver  loop  to  perturbation  of  eddy  currents  by  a  small  void  embedded  in  a 
half-space  of  uniform  electrical  conductivity  a.  Electric  current  oscillating  in  a  source 
coil,  coplanar  and  concentric  with  the  receiver  loop,  drives  eddy  currents  in  the  conducting 
half  space.  Figure  3.1  shows  the  loops  at  height  h  above  a  half  space.  A  flaw  is  at  depth 
d  and  radial  distance  rf  from  the  center  of  the  loops. 

We  first  express  the  distribution  of  eddy  currents  excited  by  the  source  in  a 
uniformly  conducting  half-space,  then  give  the  current  perturbation  produced  by  a  flaw, 
and,  finally,  determine  the  response  of  the  receiver  loop  to  perturbed  eddy  currents. 

3.1  INDUCED  EDDY  CURRENTS 

Faraday’s  law  in  the  form 


V  *  E  =  -flB/flt  (3.1a) 

-4  -♦ 

gives  the  electric  field,  E,  induced  by  a  changing  density  of  magnetic  flux,  B,  from  an 
oscillating  current.  Ampere’s  law,  expressed  by  the  relation, 

V  *  H  =  J  +  Js  (3.1b) 

4  4 

relates  magnetic  field  intensity,  H,  to  electric  current  density,  J,  in  the  conducting  region 
and  current  density,  Js  ,  of  a  source.  Equations  3.1a  and  3.1b  together  with  Ohm’s  law, 

4  4  4  ^ 

J  =  a  E,  which  relates  current  density  to  electric  field,  and  the  relation  B  =  /Xq  H  ,  with 

Ho  =  4-k  *  10 '7  H/m  ,  which  says  that  density  of  magnetic  flux  is  proportional  to  magnetic 
field  intensity,  give  current  density  and  electric  and  magnetic  fields  above  and  within  the 
conductor.  Electric  and  magnetic  fields  parallel  to  the  conductor  boundary  are  continuous 
across  the  boundary. 

To  determine  the  fields,  we  express  magnetic  flux  density  as  the  curl  of  a  vector 

potential,  A,  so  B  =  V  *  A  with  V*A  =  0.  Equations  3.1a  and  3.1b  then  tell  us  that 

V2  A  -Hoa  dA/dt  =  -Ho  Js  (3.2a) 


with  E  =  -dA/dt  and  J  =  -a  dA/dt. 

Appendix  A  derives  the  vector  potential  of  an  electric  current  oscillating  at 
frequency  w  in  a  circular  loop  of  radius  r0  at  height  h  above  a  conducting  half-space.  It 
shows  that  the  relation 


oo 

A(r,z)  =  A„(r,z)  +  (j>  f  [j=^]f]  e~k(z+h)  J  ,(kr0)  J  ,(kr)  dk  ,  (3.2b) 

gives  the  vector  potential  for  z  >  0,  where 
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z 


Figure  3.1 

Coplanar  source  coil  and  receiver  loop  at  height  h  above  a  half— space  of  uniform 
conductivtiy  a.  A  small  void  is  at  depth  d  and  radius  rf  from  the  center  of  the 
concentric  loops.  The  origin  of  cylindrical  coordinates  r,  (j>,  z  is  on  the  common  axis  of  the 
loops  at  the  surface  of  the  half-space. 
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AoM  =  [( 1-p >/2)  K (p)  -  E(p)]  |  (3.2c) 

with  p1  =  4rr0/[(z-h)J+(ro+r)2J  ,  is  the  vector  potential  of  a  current  loop  in  free  space 
(Smythe,  1969).  The  relation 

-»  „  Fr  07z~kh  i 

A(r,z)  =  polio  |  J  [  +  ij  J i(kr0)  Ji(kr)  dk  (3.2d) 

gives  it  for  z  <  0.  Here,  ?3  =  k3  +  ipQau  and  J  i(x)  is  a  Bessel  function  of  the  first  kind 

and  order  one.  Because  the  source  current  is  in  the  direction  $  of  cylindrical  coordinates 
r,  (ji,  and  z,  so  is  its  vector  potential  for  a  uniformly  conducting  naif-space. 

Electric  field  and  current  distributions  are  proportional  to  A(r,z);  namely, 

-*  ■*  4  4 

E(r,z)  =  -iu  A(r,z)  and  J(r,z)  =  -i ua  A(r,z)  ,  where  Equations  3.2b  and  3.2d  express 

4  4 

A(r,z).  Appendix  A  gives  expressions  for  the  corresponding  density  of  magnetic  flux,  B. 

3.2  EDDY  CURRENT  PERTURBATION 

Eddy  currents  induced  in  a  uniformly  conducting  half-space  by  a  current  oscillating 
in  a  circular  loop  above  the  surface  flow  in  concentric  circles  about  the  axis  of  the  source 
loop.  Following  Burrows,  1964,  we  take  a  flaw  as  a  small,  spherical  void  that  perturbs  the 
circular  flow  of  eddy  currents.  For  a  flaw  that  is  much  smaller  than  a  skin  depth  and  that 
is  far  from  boundaries,  induced  current  is  effectively  uniform  over  its  dimensions. 

A  spherical  void  imbedded  in  a  uniform  current  flow  diverts  current  around  the 
sphere.  The  diverted  flow  is  the  sum  of  a  uniform  flow  and  flow  from  a  current  dipole  at 
the  center  of  the  spherical  void.  The  relation  (Panofsky  and  Phillips,  1962) 

A  =  -  (3/2)V  J0  (3.3a) 

expresses  the  current  dipole,  A  ,  in  terms  of  the  incident,  uniform  current  J0  and  volume 
of  the  void,  V.  The  flow  outside  a  flaw  is  the  sum  J0+Ji  ,  where  the  expression 

-  f  «  (fll  -  f)]  (3.3b) 

gives  the  dipole  form  of  current  scattered  by  a  flaw,  with  ?  a  radial  vector  measured  from 
the  center  of  the  flaw. 

Current  inside  a  flaw  vanishes,  but  the  electric  field  within  a  flaw  is  (3/2)E0  > 
where  E0  =  J0/ o  and  is  the  uniform  electric  field  induced  at  a  flaw.  Total  electric  field 
outside  a  flaw  is  the  sum  E0+Ej ,  where  the  electric  field  scattered  by  a  flaw,  Ej ,  is  Ji/a . 
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The  dipole  opposes  the  induced  current  in  order  to  null  electric  current  within  a 

flaw.  It  arises  from  a  uniform  current  density  —  J0  throughout  the  sphere  that  cancels  the 
uniform  current  induced  within  the  flaw.  For  either  an  irregular  flaw  volume  or  a 
nonuniform  incident  current,  the  diverted  or  scattered  current  outside  a  flaw  comprises 
multipoles  of  higher  order. 

3.3  MAGNETIC  FLUX  PERTURBATION 

To  determine  the  perturbation  in  magnetic  flux  resulting  from  electric  current 
diverted  by  a  flaw,  w  j  again  follow  Burrows,  1964,  and  use  a  reciprocity  theorem 
(Appendix  C).  The  expression 


(3.4) 


states  reciprocity  in  terms  of  vector  potentials  and  impressed  current  densities.  It  says 
that  the  volume  integral  of  the  scalar  product  of  a  current  density  in  a  receiver  loop  ,  Jr  , 
with  the  vector  potential  from  the  current  dipole  of  a  flaw  ,  Af  ,  is  equivalent  to  the 
volume  integral  of  the  scalar  product  of  the  current  dipole  of  a  flaw,  Jf  ,  with  the  vector 
potential  from  current  in  the  receiver  loop,  Ar . 

Because  the  receiver  loop  is  a  current  filament,  the  volume  integral  containing  the 
receiver  current  density  is  simply  the  product  of  current  in  the  loop,  Ir  ,  and  flux 
threading  the  loop,  if  ,  that  results  from  the  current  dipole  of  the  flaw;  namely, 


/ 


Af.Jrd?  =  Ir#f  . 


(3.5a) 


The  vector  potential  within  a  flaw  from  current  in  the  receiver  loop,  Ar  ,  is  (3/2)A?  , 

where  A?  is  the  vector  potential  induced  at  the  flaw  from  current  in  the  receiver  loop. 

Equation  3.3a  gives  current  density  in  a  flaw,  Jf  =  -J0  ,  in  terms  of  the  current  dipole,  so 
the  relation 


•  jf  d*  =  A?- A 


(3.5b) 


expresses  the  volume  integral  in  terms  of  the  current  dipole  moment  and  the  vector 
potential  induced  at  the  flaw  by  a  current  Ir  in  the  receiver  loop.  Equations  3.4,  3.5a,  and 
3.5b  then  tell  us  that  the  expression 


*f  = 


A- A? 

~rr 


(3.5c) 


jves  the  perturbed  flux  threading  a  receiver  loop  as  a  result  of  eddy  currents  diverted  by  a 
law. 
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3.4  ELEMENTAL  FLAW  SIGNAL 


To  determine  the  signal  in  a  single  receiver  loop  from  a  flaw  illuminated  by  a  single 
source  coil,  we  combine  Equation  3.3a,  giving  the  current  dipole  that  describes  diversion  of 
eddy  currents,  with  Equation  3.5c,  giving  the  flux  threading  a  receiver  loop  from  eddy 
currents  diverted  by  a  flaw.  Namely,  the  expression 


A  =  (3/2)V  iua  Ag  ,  (3.6a) 

gives  the  current  dipole  in  terms  of  the  vector  potential  As  produced  at  a  flaw  by  a  source 
current  Is  oscillating  at  frequency  u  ,  so  Equation  3.5c  tells  us  that  the  expression 


M? 

Trrr. 


(3.6b) 


gives  flux  threading  the  receiver  loop,  #r  ,  for  unit  current  in  the  source  coil,  where  A? 
is  the  vector  potential  produced  at  the  flaw  by  a  current  Ir  in  the  receiver  loop.  Vector 
potentials  take  the  form  expressed  by  Equation  3. 2d. 


For  a  coplanar  and  concentric  source  coil  and  receiver  loop  at  height  h  above  a 
conducting  half-space,  we  use  Equation  3.2d  to  express  Equation  3.6b  in  the  form 


F(rf,d)  =  Ho  j  V  \Hq ua  Nor0  Ntri  Y0(rf,d)  Yi(rf,d)  (3.7a) 

that  gives  the  flux  threading  the  receiver  loop,  F(rf,d)  ,  resulting  from  a  flaw  of  volume  V 
at  depth  d  and  radial  distance  rf  from  the  center  of  the  loops,  for  a  unit  current  in  the 
source  coil.  The  source  coil  comprises  No  turns  of  radius  r0  ,  and  the  receiver  has  Ni 
turns  of  radius  n  .  Here,  the  expression 


f  r  -7d-kh  , 

Yi(r,d)  =  J  [  */k  +  1  j  J'(kri)  dk  <3-7b) 

gives  the  distribution  of  eddy  currents  with  depth,  d  ,  and  radial  distance,  r  ,  from  a 
driving  loop  of  radius  ri  at  height  h  .  It  measures  sensitivity  of  the  loops  to  a  flaw  at  depth 
d  and  radial  distance  r  . 
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4.0  MICROPROBE  RESPONSE  TO  A  FLAW 


Equation  3.7a  gives  the  elements  of  the  signal  from  a  flaw.  It  describes  the  response 
of  a  single  receiver  loop  to  a  flaw  illuminated  by  eddy  currents  driven  by  a  single  source 
coil.  A  microprobe  comprises  a  source  formed  by  two  opposed  coils  and  a  receiver  of  n 
loops  forming  a  gradiometer  of  order  n  .  To  describe  its  response  to  a  flaw,  we  write  the 
source  and  receiver  factors  of  Equation  3.7a  as  sums  of  factors  for  multiple  source  coils  and 
receiver  loops  of  a  microprobe.  We  number  the  two  source  coils  0  and  1  and  number  the 
n  receiver  loops  from  2  through  2+n. 

As  Figure  1.1  shows,  two  coplanar  source  coils  encircle  n  coplanar  and  concentric 
loops  of  the  receiver  that  form  a  gradiometer  of  order  n  .  Electric  current  Ii  oscillating  in 
the  inner  source  coil  balances  current  Io  oscillating  in  the  outer  coil  in  order  to  null 
interference  from  the  source  at  the  receiver.  A  balance  factor  0  gives  their  complex  ratio; 
namely,  Ij  =  —0lo  . 

The  expression 

As(rf,d)  =  /xolo  [N0r0  Y0(rf,d)  -  0N,r,  Y,(rf,d)]  (4.1) 

then  gives  the  vector  potential  of  the  source  at  a  flaw  located  at  depth  d  and  radial 
position  rf  .  Here,  Equation  3.7b  gives  the  distribution,  Yi(r,d)  ,  of  eddy  currents  with 

depth,  d  ,  and  radial  distance,  r  ,  from  a  driving  loop  of  radius  ri  at  height  h  ,  and 
Equation  2.6a  expresses  the  balance  factor  0  in  terms  of  mutual  inductances  between 
source  coils  and  receiver  loops.  Amplitude  and  phase  of  the  factor  0  specify  current  in  the 
inner  drive  coil  needed  to  make  net  flux  from  the  source  vanish  in  the  receiver  loops. 

The  expression 


■»  n  ♦  2 

[4f]n  =  M£N,riYi(r(,d)  (4-2) 

i  *2 

gives  the  net  vector  potential  at  a  flaw  for  unit  currents  in  the  n  receiver  loops.  It  tells 
sensitivity  of  a  receiver  of  order  n  to  flux  perturbed  by  a  flaw.  Equation  2.2b  specifies 
number  of  turns  Ni  and  radii  ri  of  the  n  receiver  loops  that  give  a  gradiometer  of  order 
n  . 

Equations  4.1  and  4.2  together  with  Equation  3.6b  then  tell  us  that  the  expression 

n  +  2 

Gn(rf,d)  =  Ho  j  V  i/XotrwIo^Noro  Y0(rf,d)  -  /?Nir,  Yi(rf,d)]  ^  Yi(rf,d)  (4.3) 

i  =  2 

gives  the  response  of  a  microprobe,  Gn(rf,d)  ,  with  a  receiver  forming  a  gradiometer  of 

order  n  ,  to  a  spherical  flaw  of  volume  V  at  depth  d  and  radial  position  rf  from  the 

center  of  the  microprobe.  The  response  is  the  net  flux  threading  the  receiver  for  a  current 
Io  in  the  outermost  coil  of  two  balanced  source  coils. 
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We  describe  the  response  for  microprobes  of  order  two  and  four.  A  gradiometer  of 
odd  order  gives  the  same  response  as  the  preceding  one  of  even  order,  because  the  response 
to  a  dipole  contains  even  powers  of  r  alone.  The  microprobe  carries  a  current  Io  in  the 
main  drive  coil,  oscillating  at  frequency  f  ,  and  sits  at  height  h  above  a  thick  aluminum 
plate,  with  a  =  3.54  »  107  mhos/m  .  A  spherical  flaw  0.1  mm  in  radius  is  at  depth  d  in 
the  plate  and  at  radial  distance  rf  from  the  center  of  the  microprobe.  In  each  case,  we 
give  the  net  flux  threading  the  receiver  in  units  of  a  flux  quantum,  4>0  .  for  a  current  of  one 
ampere  in  the  source;  namely,  Tn  =  Gn/<j>o/Io  >  where  Equation  4.3  gives  the  net  flux,  Gn  , 
threading  the  receiver.  A  flux  quantum  is  2.07  *  10 "l5  Tm2 . 

4.1  RESPONSE  OF  A  SECOND  ORDER  MICROPROBE 

We  give  the  response  of  a  microprobe,  with  a  receiver  formed  by  a  gradiometer  of 
second  order,  to  a  spherical  flaw  0.1  mm  in  radius  in  an  aluminum  plate  using  values  listed 
in  Table  4.1  for  number  of  turns  and  radii,  in  millimeters,  of  source  coils  and  receiver 
pickup  loops.  Source  coils  are  numbers  0  and  1 ;  receiver  loops  are  numbers  2  and  3  . 

TABLE  4.1 

COIL  SPECIFICATIONS  FOR  A  MICROPROBE  OF  SECOND  ORDER 

i  0  12  3 

Nj  16—4  1—4 

,  mm  3  2  1  0.5 


Figure  4.1  shows  magnitude  of  the  complex  response  of  a  microprobe  of  second 
order,  as  a  function  of  radial  distance,  to  a  0.1  mm  flaw  at  fixed  depths  spaced  0.5  mm 
apart,  ranging  from  0.5  mm  to  5  mm.  The  microprobe  carries  a  current  of  one  ampere  in 
the  main  drive  coil,  oscillating  at  a  frequency  of  1  kHz,  and  sits  1  mm  above  the  surface. 
Skin  depth  in  aluminum  is  2.68  mm  at  a  frequency  of  1  kHz. 

The  response  shows  a  peak  near  a  radius  of  1  mm  that  migrates  to  larger  radii  as 
flaw  depth  increases.  Peak  magnitude  decreases  markedly  at  radial  distances  and/or  flaw 

depths  greater  than  a  skin  depth,  but  still  exceeds  SQUID  sensitivity  of  10**  <f>0/4~H'z  out  to 
about  13  mm  in  radius  and  8  mm  in  depth,  as  shown  in  the  next  figure. 

Figure  4.2  shows  contours  of  constant  magnitude  in  the  r,z  plane  for  flaw  depths  to 
8  mm  and  radii  out  to  14  mm.  Contours  range  from  a  magnitude  of  10  <J>0/A  to  10’4  <j>0/A  . 
Contours  show  that  peak  response  exceeds  SQUID  sensitivity  out  to  radii  of  13  mm  or  so 
and  depths  to  8  mm,  for  a  drive  current  of  1  A.  The  line  cutting  across  contours  traces 
migration  of  peak  response  with  flaw  depth.  It  gives  a  measure  of  horizontal  resolution  of  a 
microprobe,  as  section  5.4  shows. 

figure  4.3  shows  magnitude  at  the  peak  as  a  function  of  drive  frequency,  for  flaw 
depths  of  1  mm,  3  mm,  5  mm,  7  mm,  and  10  mm.  It  shows  that  peak  response  increases 
linearly  with  drive  frequency,  reaches  a  broad  maximum,  and  then  decreases  sharply. 
Because  amplitudes  of  eddy  currents  increase  in  proportion  to  frequency,  response  increases 
linearly  with  frequency  until  losses  in  the  conductor  dominate  and  then  decrease  sharply,  as 
skin  depth  decreases  with  increasing  frequency.  For  internal  flaws  at  depths  exceeding 
1  mm,  operating  at  a  several  hundred  to  a  few  thousand  Hertz  is  best.  For  flaws  near  the 
surface,  frequencies  of  several  thousand  Hertz  are  best. 
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Figure  4.1 

Magnitude  of  the  complex  response  of  a  microprobe  of  second  order,  shown  as  a 
function  ot  radial  distance,  to  a  0.1  mm  flaw  at  fixed  depths  spaced  0.5  mm  apart,  ranging 
from  0.5  mm  to  5  mm.  The  microprobe  carries  a  current  of  one  ampere  in  the  main  drive 
coil,  oscillating  at  a  frequency  of  1  kHz,  and  sits  1  mm  above  the  surface.  Abscissa 
measures  radial  distance  in  millimeters,  and  the  left  hand  ordinate  measures  depth  in 
millimeters.  Right  hand  ordinate  measures  gain  of  the  response  in  $0/A  .  Peak  response 
for  a  0.1  mm  flaw  at  a  depth  of  0.5  mm  is  4.84  $o/A. 


r 


Figure  4.2 

Contours  of  constant  magnitude  in  the  r,z  plane  for  flaw  depths,  z,  from  zero  to  8 
mm  and  radii,  r,  from  zero  to  14  mm.  Contours  range  from  a  magnitude  of  10  $0/ A  to 
10-4  $o/A  ,  spaced  logarithmically  at  intervals  of  0.1  in  the  exponent.  The  line  cutting 
across  contours  traces  the  locus  of  peak  magnitude.  Tick  marks  on  the  locus  mark  contours 
at  decades  from  10  $0/A  to  10'4  $0/A  . 
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Figure  4.3 

Peak  response  of  a  microprobe  of  second  order  shown  as  a  function  of  frequency  for  a 
spherical  flaw  0.1  mm  in  radius  at  depths,  d  ,  of  1  mm,  3  mm,  5  mm,  7  mm,  and  10  mm  in 
a  thick  aluminum  plate.  Height  of  the  microprobe  above  the  surface,  h  ,  is  1  mm.  The 
ordinate  axis  measures  response  in  flux  quanta,  $0  ,  for  a  current  of  one  ampere  oscillating 
in  the  outermost  drive  coil  at  frequencies  ranging  from  ten  Hertz  to  ten  thousand  Hertz,  as 
marked  on  the  abscissa. 
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4.2  RESPONSE  OF  A  FOURTH  ORDER  MICROPROBE 


We  describe  response  of  a  microprobe,  with  a  receiver  formed  by  a  gradiometer  of 
fourth  order,  to  a  spherical  flaw  0.1  mm  in  radius  in  an  aluminum  plate  using  values  listed 
in  Table  4.2  for  number  of  turns  and  radii,  in  millimeters,  of  source  coils  and  receiver 
pickup  loops.  Source  coils  are  numbers  0  and  1  ;  receiver  loops  are  numbers  2,  3,  4,  and 
5  . 


TABLE  4.2 

COIL  SPECIFICATIONS  FOR  A  MICROPROBE  OF  FOURTH  ORDER 

i  0  1  2  3  4  5 

Nj  16—4  1—2  2  -2 

r* ,  mm  3  2  1  0.90  0.64  0.31 

Figure  4.4  shows  magnitude  of  the  complex  response  of  a  microprobe  of  fourth  order, 
as  a  function  of  radial  distance,  to  a  0.1  mm  flaw  at  fixed  depths  spaced  0.5  mm  apart, 
ranging  from  0.5  mm  to  3  mm.  The  microprobe  again  carries  a  current  of  one  ampere  in 
the  main  drive  coil,  oscillating  at  a  frequency  of  1  kHz,  and  sits  1  mm  above  the  surface. 

The  response  shows  a  primary  peak  near  a  radius  of  0.6  mm  and  a  secondary  peak 
near  a  radius  of  1.8  mm.  Peaks  migrate  to  larger  radii  as  flaw  depth  increase..  The  two 
peaks  form  two  lobes  on  contours  shown  in  the  following  figure. 

Figure  4.5  shows  contours  of  constant  magnitude  in  the  r,z  plane  for  flaw  depths  to 
4  mm  and  radii  out  to  5  mm.  Contours  range  from  a  magnitude  of  0.5  <J>0/A  to  10*<  <j>0/A  . 
Contours  show  that  peak  response  exceeds  SQUID  sensitivity  out  to  radii  of  about  5  mm 
and  depths  to  3.5  mm,  for  a  drive  current  of  1  A.  The  line  cutting  across  contours  traces 
migration  of  peak  response  with  flaw  depth.  It  gives  a  measure  of  horizontal  resolution  of  a 
microprobe,  as  section  5.4  shows. 


Figure  4.6  shows  magnitude  of  the  primary  peak  as  a  function  of  drive  frequency,  for 
flaw  depths  of  0.5  mm,  1  mm,  2  mm,  and  3  mm.  Peak  response  again  increases  linearly 
with  drive  frequency,  reaches  a  broad  maximum,  and  then  decreases  sharply  as  losses  in  the 
conductor  dominate.  For  internal  flaws  at  depths  exceeding  1  mm,  operating  at  a  few 
hundred  to  a  thousand  Hertz  is  best.  For  flaws  near  the  surface,  frequencies  of  several 
thousand  Hertz  are  best. 
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Figure  4.4 

Magnitude  of  the  complex  response  of  a  microprobe  of  fourth  order,  shown  as  a 
function  ofradial  distance,  to  a  0.1  mm  flaw  at  fixed  depths  spaced  0.3  mm  apart,  ranging 
from  0.5  mm  to  3  mm.  The  microprobe  carries  a  current  of  one  ampere  in  the  main  drive 
coil,  oscillating  at  a  frequency  of  1  kHz,  and  sits  1  mm  above  the  surface.  Abscissa 
measures  radial  distance  in  millimeters,  and  the  left  hand  ordinate  measures  depth  in 
millimeters.  Right  hand  ordinate  shows  gain  of  the  response,  which  is  0.025  <f>0/A  per 
division.  Peak  response  for  a  0.1  mm  flaw  at  a  depth  of  0.5  mm  is  0.108  $«,/ A. 
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Figure  4.5 

Contours  of  constant  magnitude  in  the  r,z  plane  for  flaw  depths,  z,  to  4  mm  and 
radii,  r,  to  5  mm.  Contours  range  from  a  magnitude  of  0.468  $o/A  to  10-4  $0/A  , 
spaced  logarithmically  at  intervals  of  0.1  in  the  exponent.  The  line  cutting  across  contours 
traces  the  locus  of  peak  magnitude.  Tick  marks  on  the  locus  mark  contours  in  decades 
from  10*1  $0/A  to  10*4  $o/A. 
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Figure  4.6 

Peak  response  of  a  microprobe  of  fourth  order  shown  as  a  function  of  frequency  for  a 
spherical  flaw  0.1  mm  in  radius  at  depths,  d  ,  of  0.5  mm,  1  mm,  2  mm,  and  3  mm  in  a 
thick  aluminum  plate.  Height  of  the  microprobe  above  the  surface,  h  ,  is  1  mm.  The 
ordinate  axis  measures  response  in  flux  quanta,  $0  ,  for  a  current  of  one  ampere  oscillating 
in  the  outermost  drive  c oil  at  frequencies  ranging  from  ten  Hertz  to  ten  thousand  Hertz,  as 
marked  on  the  abscissa. 
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5.0  ELECTROMAGNETIC  MICROSCOPE 


High  sensitivity  of  a  superconductive  microprobe  provides  high  resolution  needed  to 
image  internal  flaws  in  aluminum.  To  form  images,  we  cluster  microprobes  to  form  an 
array  much  like  rods  in  the  retina  of  a  magnetic  eye.  The  eye  enables  high  scan  rates  at 
high  resolution  needed  to  develop  an  electromagnetic  microscope  for  nondestructive 
evaluation  of  materials. 

Here,  we  set  forth  the  concept  of  an  array  of  microprobes  forming  a  magnetic  eye 
and  examine  its  depth  of  field,  field  of  view,  and  resolution. 

5.1  ELEMENTS  OF  A  MAGNETIC  EYE 

An  array  of  microprobes  forming  the  eye  of  an  electromagnetic  microscope 
comprises  pairs  of  parallel  rows.  Figure  5.1  shows  a  segment  of  one  pair  comprising  biaxial 
microprobes  of  second  order.  Spacing  between  centers  of  microprobes  in  each  row  is  four 
times  the  radial  field  of  view.  Offset  and  spacing  between  rows  is  set  so  microprobes  form 
a  pattern  of  equilateral  triangles,  with  their  centers  separated  by  four  times  the  radial  field 
of  view.  It  is  both  the  smallest  spacing  without  overlap  and  the  widest  spacing  that  gives 
complete  coverage. 

Closest  packing  of  microprobes  at  a  separation  of  twice  the  field  of  view  gives  double 
coverage  for  two  parallel  rows  or  complete  coverage  for  a  single  row.  Nonetheless,  closest 
packing  gives  highest  interference  between  probes.  Effectiveness  of  the  interior  source  coils 
of  each  probe  in  nulling  interference  determines  minimum  allowed  array  spacing.  Because 
of  symmetry,  null  settings  of  currents  are  equal  in  each  microprobe,  so  they  connect  in 
series. 


Because  field  of  view  changes  with  flaw  depth  and  operating  frequency,  spacing  of  a 
pair  of  rows,  or  a  single  row,  is  set  for  greatest  sensitivity  at  a  specified  depth  range  and 
operating  frequency.  A  number  of  pairs  comprises  an  eye,  each  with  a  spacing  set  for  flaws 
at  a  specified  depth.  Pairs  set  to  scan  shallow  flaws  use  a  fourth  order  gradiometer  in  order 
to  obtain  high  resolution.  High  resolution  of  a  gradiometer  of  fourth  order  limits  depth  of 
field,  so  pairs  set  to  scan  flaws  at  depths  of  a  several  millimeters  use  gradiometers  of  second 
order.  The  receiver  of  each  microprobe  forms  a  biaxial  gradiometer  in  order  to  give 
angular  resolution. 

5.2  DEPTH  OF  FIELD 

We  take  depth  of  field  as  the  depth  at  which  the  peak  signal  from  a  flaw  falls  below 

10'4  <J»o /JHz  ,  which  is  the  sensitivity  of  a  SQUID  sensor.  Depth  of  field  depends  on  flaw 
size,  gradiometer  order,  and  frequency  and  amplitude  of  the  illuminating  source  current. 

Figures  4.2  and  4.5  show  contours  of  signal  amplitude  for  gradiometers  of  second 
and  fourth  order,  respectively.  Contours  tell  depth  of  field  to  a  flaw  0.1  mm  in  radius  at 
an  operating  frequency  of  1  kHz  with  a  source  current  of  1  A.  They  show  that  depth  of 
field  is  about  8  mm  ,  for  a  gradiometer  of  second  order,  and  is  about  3.5  mm  ,  for  a 
gradiometer  of  fourth  order.  A  second  order  gradiometer  gives  somewhat  more  than  a 
twofold  advantage  in  depth  of  field. 
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Figure  5.1 

Segment  of  a  pair  of  parallel  rows  of  an  array  of  biaxial  microprobes  forming  a 
section  of  the  eye  of  an  electromagnetic  microscope.  Dashed  circles  mark  the  radial  field  of 
view  for  each  microprobe.  Dashed  lines  connecting  their  centers  highlight  the  equilateral, 
triangular  pattern,  with  a  spacing  between  centers  of  four  times  the  radial  field  of  view. 
The  arrow  marks  the  direction  of  scan,  and  dots  denote  continuing  replication  of  the  array. 
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Figures  4.3  and  4.6  give  peak  signal  amplitude  as  a  function  of  operating  frequency 
for  flaws  at  fixed  depths,  tor  gradiometers  of  second  and  fourth  order  respectively.  They 
show  that  drive  frequencies  of  a  few  kilohertz  give  greatest  signal  amplitudes  for  flaw 
depths  of  a  few  millimeters.  Depth  of  field,  then,  is  greatest  at  frequencies  of  a  few 
kilohertz  and  changes  little  with  frequency  at  frequencies  around  1  kHz. 

Because  signal  amplitude  is  proportional  to  volume  of  a  flaw,  depth  of  field  increases 
as  the  cube  of  flaw  size.  It  also  increases  linearly  with  amplitude  of  the  drive  current, 
provided  interference  from  the  source  remains  null. 

5.3  FIELD  OF  VIEW 

We  take  field  of  view  as  the  maximum  radius  at  which  the  signal  from  a  flaw  falls 

below  10 '4  ♦o/JHi  ,  for  a  flaw  at  a  fixed  depth.  It  depends  on  flaw  size,  gradiometer  order, 
and  frequency  and  amplitude  of  the  drive  current,  as  does  depth  of  field. 

Contours  at  10  "4  $0/A  in  Figures  4.2  and  4.5  trace  field  of  view  as  a  function  of  flaw 
depth,  for  second  and  fourth  order  gradiometers,  respectively,  operating  at  a  frequency  of 
1  kHz  with  a  drive  current  of  1  A.  Flaw  size  is  again  0.1  mm.  They  show  that  field  of 
view  varies  widely  with  depth  of  a  flaw. 

For  a  gradiometer  of  second  order,  field  of  view  is  about  3  mm  for  a  flaw  at  the 
depth  of  field  of  about  8  mm.  It  increases  to  nearly  13  m  for  flaws  near  the  surface.  For  a 
gradiometer  of  fourth  order,  field  of  view  is  about  1.5  mm  for  a  flaw  at  the  depth  of  field  of 
about  3.4  mm.  At  greater  depths,  field  of  view  vanishes.  At  shallower  depths,  it  increases 
over  the  first  lobe  of  the  contour  to  about  3  mm  at  a  depth  of  about  2.4  mm.  Over  the 
second  lobe,  it  increases  to  a  maximum  value  of  about  5  mm  at  a  depth  of  about  1  mm. 

5.4  RESOLUTION 

Contours  of  constant  signal  amplitude  in  a  horizontal  plane  above  a  flaw  are  circles 
centered  on  the  flaw.  They  give  an  image  of  a  flaw  that  results  from  repeated  scans  over  a 
flaw  of  a  single  microprobe  or  a  single  scan  of  a  magnetic  eye. 

Figure  5.2a  and  Figure  5.2b  show  them  at  a  height  of  1  mm  for  a  0.1  mm  flaw  at  a 
depth  of  1.44  mm  and  4.96  mm,  respectively,  for  a  gradiometer  of  second  order  with  a  drive 
current  of  1  A  oscillating  at  1  kHz. 

figure  5.3a  and  Figure  5.3b  show  them  for  a  gradiometer  of  fourth  order  at  depths 
of  0.96  mm  and  3.04  mm,  respectively,  again  at  a  height  of  1  mm  for  a  0.1  mm  flaw  with  a 
drive  current  of  1  A  oscillating  at  1  kHz. 

An  eye  gives  multiple  images,  one  for  each  row  or  pair  of  rows.  They  depend  on 
gradiometer  order  and  microprobe  spacing  in  a  row  as  well  as  frequency  and  amplitude  of 
the  drive  current  of  a  row.  Multiple  images  can  localize  depth  of  a  flaw.  Gradiometers  of 
fourth  or  higher  order  identify  shallow  flaws,  and  gradiometers  of  second  order  see  deep  as 
well  as  shallow  flaws.  For  shallow  flaws,  ratios  of  image  amplitudes  from  gradiometers  of 
second  and  fourth  order,  for  example,  are  independent  of  flaw  size  and  depend  on  flaw 
depth  alone,  for  fixed  frequency  and  amplitude  of  drive  currents.  The  information  gives 
means  of  locating  depth  of  a  flaw.  Alternatively,  changing  drive  current  and  or  frequency 
of  a  row  moves  its  depth  of  field  to  the  of  depth  of  a  flaw  and  so  determines  its  depth. 
Using  a  pulsed  waveform  can  also  give  flaw  depth  from  the  time  delay  of  the  scattered 
field. 
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Figure  5.2b 

Contours  of  constant  signal  amplitude  in  a  horizontal  plane  at  a  height  of  1  mm  for 
a  0.1  mm  flaw  at  a  depth  of  4-96  mm  for  a  gradiometer  of  second  order  operating  with  a 
drive  current  of  1  A  oscillating  at  1  kHz.  The  dashed  circle  marks  the  peak  signal 
amplitude.  It  gives  a  measure  ofhorizontal  resolution. 
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Contours  of  constant  signal  amplitude  in 
a  0.1  mm  flaw  at  a  depth  of  3.04  mm  for  a  gra 
drive  current  of  1  A  oscillating  at  1  kHz  ' 
amplitude.  It  gives  a  measure  of  horizontal 


Radius  of  the  contour  marking  the  peak  signal,  highlighted  by  dashed  circles  in 
Figures  5.2a,  5.2b,  5.3a,  and  5.3b,  gives  a  measure  of  horizontal  resolution  for  a  magnetic 
eye.  Radius  increases,  so  resolution  decreases,  linearly  with  depth  of  a  flaw.  Figures  5.2a 
and  5.2b  show  the  decrease  in  resolution  from  a  flaw  depth  of  1.44  mm  to  4.96  mm,  for  a 
gradiometer  of  second  order.  Figures  5.3a  and  5.3b  show  it  from  0.96  mm  to  3.04  mm,  for 
a  gradiometer  of  fourth  order. 

The  locus  of  peak  signal  in  Figures  4.2  and  4.5  gives  the  linear  change  in  resolution 
with  flaw  depth,  for  gradiometers  of  second  and  fourth  orders  respectively.  For  second 
order,  we  find  that  the  relation 


rp  =  0.371  d  4-  0.792  ,  millimeters,  (5.1a) 

gives  the  radius  of  the  contour  marking  the  peak  signal,  rp  ,  for  a  flaw  depth  d  ,  both 
measured  in  millimeters.  Radius  increases  linearly  from  0.792  mm  for  a  surface  flaw  to 
3.76  mm  at  the  depth  of  field  8  mm.  For  fourth  order,  the  relation 

rp  =  0.256  d  +  0.65  ,  millimeters,  (5.1b) 

gives  the  contour  radius.  It  increases  linearly  from  0.65  mm  for  a  surface  flaw  to  1.55  mm 
at  the  depth  of  field  of  3.52  mm. 

Spacing  of  contours  Figures  4.2  and  4.5  give  a  measure  of  vertical  resolution  for 
gradiometers  of  second  and  fourth  orders  respectively.  Differences  between  contours 
represent  a  26  %  change  in  drive  current  amplitude.  Near  the  depth  of  field,  a  26  %  change 
in  current  amplitude  gives  a  change  in  depth  of  field  of  0.24  mm  for  a  second  order 
gradiometer  and  0.133  mm  for  a  fourth  order  gradiometer.  We  expect,  then,  a  vertical 
resolution  of  a  few  tenths  of  a  millimeter  for  a  0.1  mm  flaw. 
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6.0  ASSESSMENT  OF  FEASIBILITY  AND  TECHNICAL  RISK 

Successful  development  of  an  electromagnetic  microscope  for  eddy  current  testing 
requires  developing  a  microprobe,  forming  a  magnetic  eye  from  an  array  of  microprobes, 
and  demonstrating  that  performance  of  the  eye  surpasses  that  of  conventional  eddy  current 
probes.  Here  we  examine  feasibility  of  the  design  of  a  microprobe  and  of  a  magnetic  eye, 
identify  main  problems  facing  development,  and  estimate  technical  risk  of  developing  an 
electromagnetic  microscope. 

6.1  FEASIBILITY 

Analyses  presented  herein  show  that  design  of  a  high  performance  microprobe  is 
feasible.  A  current  of  1  A  oscillating  at  a  few  kilohertz  in  a  drive  coil  3  mm  in  radius 
excites  large  enough  eddy  currents  in  an  aluminum  plate  to  give  an  appreciable  signal  in 
pickup  loops  1  mm  in  radius  from  a  0.1  mm  flaw  at  a  depth  of  a  few  millimeters.  Peak 
signal  amplitude  is  several  flux  quanta,  for  a  second  order  microprobe,  and  a  few  tenths  of 
a  flux  quantum,  for  a  fourth  order  microprobe.  Peak  amplitudes  far  exceed  the  probe 

sensitivity  of  10*4  <j>o/JITz  .  Loop  sizes  are  well  within  the  art  of  miniature  fabrication 
techniques,  and  current  densities  required  are  modest  for  superconductive  wire. 

Although  analyses  demonstrate  feasibility  of  the  configuration  and  loop  sizes  of  a 
microprobe,  they  allow  performance  to  increase  indefinitely  with  drive  current.  In  practice, 
interference  from  the  source  coils  sets  the  maximum  allowable  drive  current  and  so  limits 
performance.  Nulling  interference  from  the  source  is  the  central  problem  facing 
development  of  a  microprobe  as  well  as  an  array  forming  a  magnetic  eye. 

Its  solution  comprises  two  main  parts.  First  is  to  control  current  in  the 
compensating  drive  coil  in  order  to  suppress  interference  in  the  pickup  loops.  Second  is  to 
apply  an  external  feedback  signal  to  the  SQUID  sensor  in  order  to  null  residual 
interference.  Interference  from  one  ampere  in  the  main  drive  coil,  3  mm  in  radius,  is  of  the 
order  of  107  flux  quanta  without  compensation.  External  feedback  to  a  SQUID  sensor  can 
compensate  residual  interference  of  a  hundred  flux  quanta  or  so.  Cuirrent  in  the 
compensating  coil,  then,  must  suppress  interference  by  a  factor  of  104  to  105  in  order  to 
reduce  it  below  100  flux  quanta.  A  circuit  using  a  phase-locked  oscillator  to  control 
amplitude  and  phase  of  current  in  the  compensating  coil  within  a  few  ppm  of  set  values  is 
within  the  art,  so  suppression  by  as  much  as  a  factor  of  10®  is  practical.  An  overall, 
combined  suppression  by  a  factor  of  10s,  then,  is  within  the  art,  so  drive  currents  of  1  to  10 
amperes  in  a  drive  coil  3  mm  in  radius  are  feasible. 

Compensation  of  a  symmetric  array  of  microprobes  follows  from  that  of  a  single 
microprobe.  Symmetry  of  an  array  says  that  compensation  for  each  microprobe  is  the 
same,  so  they  connect  in  series.  Drive  coil  compensation  for  a  symmetric  array,  then,  is  the 
same  as  for  a  single  microprobe.  Feedback  to  individual  SQUID  sensors  compensates 
residual  asymmetry  of  an  array. 

The  remaining  principal  problem  facing  development  of  a  microprobe  is  its 
cryogenic  structure.  Design  of  the  cryogenic  structure  poses  two  main  problems.  The  first 
is  to  achieve  a  spacing  of  1  mm  or  so  between  the  pickup  coils  and  a  test  surface.  The 
second  is  to  attach  the  probe  with  a  cryogenic  umbilical  to  a  remote  cryostat  housing 
SQUID  sensors.  The  first  requirement  is  essential  for  high  performance;  the  second  makes 
it  convenient  to  scan  test  samples. 
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Radiant  heat  tranfer  from  a  polished  surface  at  300  K,  with  an  emissivity  of  0.03,  to 
one  at  4  K  is  about  0.007  mW/(mm)!,  so  the  radiant  heat  load  on  the  tip  of  a  microprobe  4 
mm  in  radius  is  about  0.35  mW.  A  copper  wire  1  mm  in  radius,  with  a  thermal 
conductivity  of  about  1  W/cm/K  at  4  K,  requires  a  temperature  difference  of  1  K  over  a 
length  of  100  cm  to  transport  a  heat  load  of  0.35  mW.  Consequently,  fine  copper  wire  in  a 
flexible,  stainless  steel  vacuum  jacket,  which  is  superinsulated,  can  provide  a  cyrogenic 
umbilical  cord  one  meter  or  so  long  that  links  the  microprobe  tip  to  a  cryostat  containing 
liquid  helium. 

At  the  tip  of  the  microprobe  is  a  thin  sapphire  window  that  allows  a  1  mm  spacing 
between  the  pickup  coils  and  a  test  surface.  A  point  contact  fixes  pickup  coils  in  the  tip  to 
prevent  relative  motion  from  thermal  contraction  during  cooling  to  4  K.  Differential 
contraction  is  taken  up  by  slack  in  the  copper  wire  that  thermally  grounds  the  tip  to  the 
helium  reservoir  in  the  cryostat.  The  cryostat  houses  the  SQUID  sensor  and  is  a 
conventional  design  that  uses  thermal  shields  cooled  by  helium  vapor.  A  probe  that  fits 
down  the  central  tube  of  the  cryostat  immerses  the  sensor  in  the  helium  reservoir.  An 
umbilical  cord  leads  out  from  the  top  of  the  cryostat  probe  and  enables  readily  scanning 
surfaces  with  a  microprobe. 

We  expect  ac  losses,  which  depend  on  wire  characteristics  and  operating  frequency, 
from  a  current  of  one  ampere  in  the  superconducting  drive  coils  to  contribute  a  few 
milliwatts  to  the  heat  load  on  the  umbilical  cord,  which  would  increase  its  temperature 
differential  by  several  degrees  Kelvin.  The  increase  can  be  enough  to  exceed  the  critical 
temperature  of  niobium,  even  though  helium  loss  is  acceptable  (  a  loss  of  one  liter  of  liquid 
per  day  requires  30  mW).  Operating  at  drive  currents  of  a  few  amperes  and  frequencies  of 
several  thousand  Hertz  could  increase  the  heat  load  tenfold.  It  would  require  tripling  the 
the  diameter  of  copper  wire  in  the  umbilical  cord  or  using  a  cryogenic  heat  pipe  (Chi  and 
Cygnarowicz,  1970)  to  obtain  a  smaller  diameter  structure  with  thermal  conductance 
tenfold  greater  than  copper. 

6.2  TECHNICAL  RISK 

Technology  of  niobium  SQUIDs  enables  development  of  a  microprobe.  Reliable 
sensors  are  commercially  available.  The  main  technical  risk  is  in  meeting  performance 
goals  and  developing  a  cryogenic  design  that  makes  an  electromagnetic  microscope  a 
versatile  research  tool.  Almost  certaimy,  a  first  probe  design  falls  short  of  anticipated 
performance  and  utility,  but  an  iterative  development  plan  can  advance  design  of  a 
microprobe  and  increase  its  performance  and  utility  in  sequential  stages.  The  plan 
minimizes  risks  in  meeting  performance  goals. 

Cryogenic,  mechanical,  and  electronics  designs  present  the  main  risk  factors  in 
development.  Mechanical  and  electronics  design  are  critical  elements  that  determine 
performance  and  so  present  the  greatest  risk.  Cryogenic  design  mainly  determines  utility 
of  a  microprobe  but  also  can  affect  its  performance. 

Mecbanical  and  electronic  design  set  performance  limits  of  a  microprobe. 
Suppressing  interference  and  noise  of  a  microprobe  requires  precise  mechanical  design  and 
machining.  Electronic  feedback  nulls  interference  that  remains  because  of  limits  to 
mechanical  precision.  It  requires  high  quality  electronics  to  null  residual  interference  and 
to  maintain  a  stable  operating  point  with  minimal  drift  and  harmonic  contamination.  To 
achieve  high  performance,  combined  mechanical  precision  and  electronic  feedback  must 
reduce  interference  by  a  factor  of  107  to  10®,  with  mechanical  precision  contributing  a 
factor  of  103  to  104.  The  required  suppression  of  interference  falls  within  practical  limits  of 
high  quality  mechanics  and  electronics  and  so  presents  a  moderate  risk  io  achieving  high 
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performance. 


Development  of  a  cryogenic  umbilical  to  thermally  ground  a  microprobe  to  a 
cryostat  housing  SQUID  sensors  enhances  utility  of  a  microprobe  and  sets  a  secondary 
limit  on  its  performance.  Use  of  copper  wire  to  transport  heat  in  an  umbilical  presents  a 
moderate  risk,  but  it  could  limit  performance  of  large  arrays  of  microprobes.  Use  of  a 
helium  filled  heat  pipe  to  transport  heat  would  remove  performance  limits  for  large  arrays, 
but  it  presents  a  higher  risk,  because  it  is  a  new  cryogenic  technique. 

Overall  technical  risk  of  developing  high  peformance  microprobes  and  using  them  to 
form  the  magnetic  eye  of  an  electromagnetic  microscope  is  moderate. 
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7.0  CONCLUSION 


High  performance  microprobes  hold  the  promise  of  forming  the  magnetic  eye  of  an 
electromagnetic  microscope  that  promises  high  resolution,  wide  field  of  view,  and  large 
depth  of  field  for  imaging  internal  flaws  in  aluminum.  A  current  of  1  A  oscillating  at  a  few 
kilohertz  in  drive  coils  a  few  millimeters  in  radius  together  with  coplanar  receiver  loops 
1  mm  or  less  in  radius  form  a  high  performance  microprobe.  A  microprobe  with  a  receiver 
forming  a  gradiometer  of  fourth  order  gives  a  field  of  view  of  about  5  mm  and  a  horizontal 
resolution  of  about  1  mm  to  a  0.1  mm  flaw  at  a  depth  of  a  few  millimeters.  A  receiver 
forming  a  gradiometer  of  second  order  gives  a  field  of  view  of  12  mm  or  so  and  a  horizontal 
resolution  of  about  3  mm  at  a  depth  of  field  of  8  mm. 

Closely  packed  microprobes  arrayed  in  parallel  rows  form  the  eye  of  a 
electromagnetic  microscope.  The  eye  gives  multiple  images  of  an  internal  flaw,  one  for 
each  row  or  pair  of  rows,  that  can  resolve  depth  of  a  0.1  mm  flaw  to  a  few  tenths  of  a 
millimeter  together  with  a  horizontal  resolution  of  1  mm  or  so.  Cryogenic  umbilicals 
connect  the  eye  to  a  remote  cryostat  housing  SQUID  sensors  to  provide  ease  of  scanning 
and  make  it  a  versatile  tool. 
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APPENDIX  A 


EDDY  CURRENTS  IN  A  UNIFORMLY  CONDUCTING  HALF-SPACE 
FROM  AN  ELECTRIC  CURRENT  OSCILLATING  IN  A  CIRCULAR  LOOP 

Appendix  A  describes  eddy  currents  induced  in  a  half-space  of  uniform  electrical 
conductivity,  <7,  by  electric  current  oscillating  in  a  circular  loop  at  height  h  above  the 
surface  of  the  half-space.  The  axis  of  the  loop  of  radius  r0  is  normal  to  the  surface. 

We  use  Ampere’s  law  together  with  Faraday’s  law  to  formulate  a  description  of 

induced  currents  and  fields  in  terms  of  a  vector  potential,  A,  of  magnetic  flux  density,  B. 
We  express  the  vector  potential  in  terms  of  two  scalar  potentials  and  obtain  solutions  for 

them  in  terms  of  Fourier  integrals  in  cylindrical  coordinates  r,  (j>,  and  z,  where  the  z  axis 
points  along  the  axis  of  the  current  loop.  Scalar  potentials  give  so-called  transverse 
electric  and  transverse  magnetic  type  fields.  Source  currents  flowing  tangential  to  the 
surface  of  the  half-space  generate  transverse  electric  type  fields,  and  source  currents 
flowing  normal  to  the  surface  generate  transverse  magnetic  type  fields. 

A.1  FORMULATION 

Ampere’s  law  in  the  form 


V  «  H  =  J  +  J8  (Al.a) 

together  with  the  constituent  relation,  B  =  foR  with  =  4x  *  10 *7  H/m  ,  and  Ohm’s 
law,  J  =  <jE,  relate  magnetic  flux  density  to  electric  field,  E,  induced  by  an  imposed, 
oscillating  current  of  density  J8.  Faraday’s  law,  expressed  by  the  relation 


V  *  E  =  —  dB/dt  ,  (Al.b) 

closes  the  formulation  by  relating  the  induced  electric  field  to  a  changing  density  of 
magnetic  flux,  B. 

4  4 

Equations  Al.a  and  Al.b  tell  us  that  a  solenoidal  vector  potential,  7-A  =  0, 
satisfies  the  relation 

V2A  -  w  dk/dt  =  -fioh  ,  (A2) 


where  B  =  V  *  A  and  E  =  —dk/dt  .  Equation  A2  together  with  continuity  of  tangential 
components  of  electric  and  magnetic  fields  across  boundaries  give  the  vector  potential  for  a 
specified  source  current. 

A.2  SCALAR  POTENTIALS 

Because  the  vector  potential  is  solenoidal,  we  express  it  as  the  curl  of  a  second 

vector  potential,  Z  ,  composed  of  two  scalar  potentials  Zt  and  Z2  ;  namely,  (Smythe, 
1969) 


A1 


Z  =  Ziz  +  V*Z2z  ,  (A3. a) 

where 

A  =  /ioV  *  Z  .  (A3.b) 

Equations  A3.a  and  A3.b  tell  us  that  the  relations 

E  =  d/dt  [  z  *  VZ,  +  z  7%  -  7(£-VZ2)  ]  (A3.c) 

and 

H  =  z  *  7(72Z2)  -  z  72Z,  +  7(z-7Z,)  (A3.d) 

give  electric  and  magnetic  fields  in  terms  of  scalar  potentials.  Equation  A3.c  says  that  the 

scalar  potential  Zi  gives  an  electric  field  that  is  transverse  to  the  z  axis,  a  so-called 
transverse  electric  type  field,  and  Equation  A3.d  says  that  the  scalar  potential  Z2  gives  a 
transverse  magnetic  type  field. 

-* 

We  find  from  Equations  A2  and  A3.b  that  the  vector  potential  Z  satisfies  the 
relation 


V2Z  — /iocr  dZ/dt  =  —  Ms  ,  (A4.a) 

where  M„  is  the  magnetization  density  of  the  source;  namely, 

J8  =  7  *  Ms  .  (A4.b) 

Because  the  source  current  is  solenoidal,  we  can  express  source  magnetization  in 
terms  of  two  scalar  components,  Mj  and  M2 ,  as 

Mg  =  Mi  z  +  V  *  M2  z  ,  (A5.a) 

so 

js  =  -  z  *  VM!  -  z  7^2  +  7(z-7M2)  .  (A5.b) 

The  component  Mi  gives  electric  currents  transverse  to  z.  A  current  along  z  requires  a 
transverse  component  of  magnetization,  expressed  by  M2 . 

Equations  A3.a  and  4.a  then  tell  us  that  the  scalar  potentials  Zt  and  Z2  satisfy  the 
relation 


V2Zi  - /ioa  3Zj/dt  =  -  Mi  ,  (A6) 

with  i  =  1,  2.  Equation  A6  tells  us  that  the  scalar  source  magnetization  Mi  excites  a 
transverse  electric  type  field  and  that  M2  excites  a  transverse  magnetic  type  field. 

Consequently,  a  source  current  transverse  to  the  z  axis  excites  a  transverse  electric  type 

field,  specified  by  the  scalar  potential  Zt  ,  and  a  current  along  the  z  axis  excites  a 
transverse  magnetic  type  field,  specified  by  the  scalar  component  Z2  alone.  Current 
flowing  in  a  circular  loop  parallel  to  the  surface  of  a  half-space,  then,  excites  a  transverse 
electric  type  field  alone. 


A2 


A.3  TRANSVERSE  ELECTRIC  TYPE  FIELD  OF  A  CIRCULAR  CURRENT 


The  relation 


J8  =  <M  £(r-r0)  6(z-h)  (A7.a) 

expresses  the  current  density  in  a  circular  loop  of  radius  r0  that  carries  a  current  I  in 

direction  $  at  height  h  above  the  surface  of  a  half-space.  Dirac  delta  functions,  6(x)  , 
express  localization  of  the  current.  Equation  A5.b  tells  us  that  the  relation 

Ms  =  -  z  I  % z-h)  (A7.b) 

gives  the  corresponding  density  of  magnetization,  where  <%(x)  is  a  Heaviside  step 
function.  Consequently,  fields  of  the  current  loop  are  a  transverse  electric  type,  so 

E  =  i uz*  VZj  (A7.c) 

and 

H  =  —  z  V2Zi  -f  V(z-VZi)  (A7.d) 

for  a  current  oscillating  at  frequency  u. 

Equation  A6  then  tells  us  that  the  expressions 

V2Z^  =  I  <^(r-r0)  6(z— h)  ,  for  z  >  0,  (A8.a) 

and 


V2Zj  —  iu(i0aZcl  =s  0  ,  for  z  <  0,  (A8.b) 

specify  a  scalar  potential  above  the  surface,  Z\  ,  and  within  the  conducting  half-space, 

Z\  .  Continuity  of  tangential  components  of  electric  and  magnetic  fields  at  z  =  0 
requires  that  Zi  and  its  derivative  dZ\fdz  be  continuous  across  the  surface;  namely,  that 
Z*  =  Z\  and  dZ*/dz  =  dZcJdz  at  z  =  0.  Equations  8. a  and  8.b  together  with  continuity 

required  at  the  boundary  determine  the  scalar  potentials  Z\  and  Z\  .  We  express  the 
potentials  in  terms  of  Fourier  integrals. 

A.3.1  Fourier  Integrals 

The  integral 


Z(?,z)  =  JJ  Z(k,z)e“*'f  dk 


(A9.a) 


-00 
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where 


OD 

Z(k,z)  =  (2ir)'2  jj  Z(?,z)  e*'f  d? 

— OD 


(A9.b) 


expresses  a  function  Z(?,z)  in  terms  of  its  Fourier  transform  Z(k,z)  over  a  radial  vector 
f .  From  Equations  A8.a  and  A8.b,  we  then  find  that  the  expressions 

-  V  Z!  =  ~lr°2^kl°)  l(z-k)  ,  for  z  >  0,  ( A9.c) 

and 


j^-7,Zt  =  0  ,  for  r  <  0, 


(A9.d) 


together  with  continuity  requirements  at  z  =  0,  Z\  =  Z\  and  dZ*/dz  =  dZcJdz  ,  specify 
Fourier  transforms  of  the  scalar  potential  Z\ ,  where  72  =  k2  +  iui  . 


We  use  Equations  A9.c  and  A9.d  to  find  that  the  expressions 

+Br^e'k(z+ll)] 


and 


and 


OD 


nQ7z-kh  -1 

l+7/fc  j  Ji(kro)  Jo(kr)  dk/k 


give  corresponding  scalar  potentials.  Here,  the  relation 


— 

1  + 

l+4/(&)2  +i 

< 

l+4/(£k)2  -  1 

P 

expresses  the  complex  ratio  7/k  in  terms  of  the  skin  depth  6  =  . 

A4 


(AlO.a) 


(AlO.b) 


„A_Ir0  Ji(kro) 

zi^“2t  k2(l+7/kJ 

give  the  transforms  Z\  ,  for  z  >  0  and  Z\  ,  z  <  0  .  The  integrals 

OD 

Z»M  =  *§»/[ e-klz-hl  +  flT$$<rk(Z+h)]  Jl0“o)  Jo(kr)  dk/k  (AlO.c) 


(AlO.d) 


(AlO.e) 


A.3.1  Vector  Potential 


From  the  relation 

A  =  /ioV*Z  =  -$/Xo  dZifdi 
and  Equations  AlO.c  and  AlO.d,  we  find  that  the  expressions 


(All.a) 

(All.b) 

(All.c) 


give  the  vector  potential  Aa  ,  for  z  >  0  and  Ac  ,  for  z  <  0  .  For  low  frequencies  and  or 

low  conductivity,  7  «  k  and  Aa  s  Ac  »  A0  ,  where  A0  is  the  vector  po*ei  tial  of  a  loop  in 
free  space;  namely, 

OD 

A0(r,z)  =  J  e~k|z_h|  J,(kr0)  J^kr)  dk  ,  (All.d) 

0 

or  (Smythe,  1969) 

Ao(r,z)  =  f  [(1— pa/2)  K(,i)  -  E(p)  ]  i  ,  (All.e) 

where  K (p)  and  E(p)  are,  respectively,  complete  elliptic  integrals  of  the  first  and  second 
kinds  and  p2  =  4rr0/[(z— h)2+(r+r0)2j.  Equations  All.b  and  All.c  expressing  the  vector 
potential  agree  with  expression  given  by  Dodd,  Deeds,  and  Luquire,  1969. 

A.3.2  Electric  Field  and  Eddy  Currents 

The  relation 

E(r,z)  =  -iw  {  c^(z)  Aa(r,z)  +  [l-^(z)]  Ac(r,z)  }  (A12.a) 

gives  the  corresponding  electric  field,  where  Jiz)  denotes  a  Heaviside  step  function,  and 
the  relation 

J(r,z)  =  -iua  Ac(r,z)  (A12.b) 

gives  eddy  currents  induced  in  the  conducting  half-space,  z  <  0. 


A5 


* 


A.3.3  Magnetic  Flux  Density 
The  relations 

00 

Br(t,z)  =  Br0(r,z)  +  e_k(z+h)  J,(kr0)  J,(kt)  k  dk  (A13.a) 

and 

00 

B«(t,z)  =  Bz0(r,z)  +  ^  5k(z+h)  ]  J,(kr0)  J„(kr)  k  dk  (A13.b) 

give  respective  radial  and  vertical  components  of  the  density  of  magnetic  flux  for  z  >  0, 
and  the  relations 


and 


Br(r,z)  =  -Hoh0 


e72- kh  J,(kr0)  Ji(kr)  k  dk 


Bz(r,z)  —  /iolr0 


e7z— kh  j^kr^  J(j^kr^  k  ^ 


(A13.c) 


(A13.d) 


give  the  components  for  z  <  0.  Here,  Br0  and  Bz0  are  respective  radial  and  vertical 
components  of  the  magnetic  flux  density  of  a  loop  in  free  space;  namely, 


Br0(r,z) 


(A13.e) 


and 


Bz0(r,z) 


Ur+r0)2+(z-h)2 


+ KM  ] 


(Al3.f) 


where  again  K(p)  and  E(p)  are,  respectively,  complete  elliptic  integrals  of  the  first  and 
second  lands  and  p2  =  4rr0/l(z-h)2+(r+r0)2]. 
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APPENDIX  B 

MUTUAL  INDUCTANCES  OF  COAXIAL  LOOPS 
ABOVE  A  CONDUCTING  HALF-SPACE 

Appendix  B  uses  the  vector  potential  derived  in  Appendix  A  to  give  an  expression 
for  the  mutual  inductance  of  two  coaxial  loops  above  a  conducting  naif-space.  Mutual 
inductance  of  a  pair  of  loops,  Mij  ,  is  the  magnetic  flux  threading  one  loop  as  a  result  of  a 
unit  electric  current  flowing  in  the  other  loop. 

For  coaxial,  circular  loops  with  radii  ri  and  rj  and  corresponding  numbers  of 
turns  Ni  and  Nj  ,  the  expression 


2t 

Mij  =  Niri  Nj  / f  Aj(ri,hi)  d<|>  (Bl) 

o 

-* 

gives  the  flux  threading  the  i th  loop  for  unit  current  flowing  in  the  j th  loop,  where  Aj  is 
the  vector  potential  of  current  in  the  jth  loop.  It  is  evaluated  at  the  radius  and  height,  ri 
and  hi  ,  of  the  ith  loop.  Integration  in  Equation  BI  extends  over  the  azimuthal 
coordinate  $  of  cylindrical  coordinates  r,  <j),  and  z  measured  from  the  axis  of  the  coaxial 
loops.  For  coaxial,  circular  loops,  the  vector  potential  is  azimuthally  symmetric  and  points 

along  the  unit  vector  $  ,  so  integration  simply  gives  a  factor  of  2  r . 

Equation  All.b  of  Appendix  A  gives  the  vector  potential  Aj(r,h)  ,  h>0  ,  of  a 

current  Ij|  flowing  in  a  circular  loop  of  radius  rj  at  height  hj  above  a  conducting 
half-space;  namely, 

OD 

Aj(r,h)  =  Aj(r,h)  +  f  e-k(h+hi)  J,(ki)  J,(kij)  dk  .  (B2.a) 

0 


Here,  A°  is  the  vector  potential  of  the  loop  in  free  space;  namely, 


00 


A°(r,h)  =  e“klh-hjl  J^kr)  J,(krj)  dk  ,  (B2.b) 


or  (Smythe,  1969) 


Aj(r.h)  =  J|i  [(l-p>/2)  K(pj)  -  E(pj)  ]  $  ,  (B2.c) 

where  K(pj)  and  E(pj)  are,  respectively,  complete  elliptic  integrals  of  the  first  and  second 
kinds  and  pj  =  4rrj/[(h— hj)J+(r+rj)2]. 


Bl 


<* 

* 

• 

Equations  B1  and  B2.a,  with  Ij  =  1  ,  then  tell  us  that  the  expression 

OD 

My  =  M?j  +  N,t,  Njtj  /  feZ©e^(lli+hj>  J,(kr,)  J.fkrj)  dk 

0 

(B3.a) 

gives  the  mutual  inductance  of  a  pair  of  coaxial  loops  of  radii  ri  and  rj  and 
corresponding  numbers  of  turns  Ni  and  Nj  at  heights  hi  and  hj  above  a  conducting 

half-space.  Here,  is  the  mutual  inductance  of  the  loops  in  free  space;  namely, 

00 

M°j  =  x/io  Niri  Njrj  J  e"klhi_hil  J,(kri)  J,(krj)  dk  , 

a 

(B3.b) 

or 

My  =  PoNjNj  [(hH«j)»  +  (Wj)1]  [(l-Pij/2)  K(ptJ)  -  E(Pij)  ]  , 

(B3.c) 

where  again  K(py)  and  E(pij)  are,  respectively,  complete  elliptic  integrals  of  the  first  and 

second  kinds  and  p\ j  =  4rirj/[(hi-hj)2+(ri+rj)2].  Equation  B3.a  says  that  My  = 
required  by  reciprocity. 

Mji  ,  as 

For  coplanar  loops  at  height  h  above  the  surface,  hi  =  hj  =  h  ,  so  Equation  B3.a 
says  that 

OD 

My(h)  =  MJJ  +  m  N,I,  Njrj  /  e-2kh  J,(krj)  J.fkrj)  dk  , 

0 

(B4.a) 

where 

OD 

M?]  =  x/io  Nin  Njtj  J  Ji(kri)  Ji(krj)  dk  , 

(B4.b) 

0 

or  (Gradshteyn  and  Ryshik,  1980) 

M?J  =  PoNiNj  (ri+rj)>  {[HPy>/2]  K(Py  -  E(^)}  , 

(B4.C) 

where  (p?j) 2  =  4rirj/(ri+rj)2 . 

B2 
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APPENDIX  C 

RECIPROCITY  RELATION 

Appendix  C  develops  a  reciprocity  relation  between  a  vector  potential  Ai  that 
results  from  a  source  current  JSi  and  a  vector  potential  Aj  that  arises  from  a  source 
current  JSj  .  We  use  the  vector  identity  (Gradshteyn  and  Ryzhik,  1980) 

V-(Ai  x  Hj)  =  Hj-V  «  Ai- Ai-V  «  Hj  (Cl.a) 

together  with  Ampere’s  law  in  the  form 

V  x  Hj  =  Jj  +  Jsj  (Cl.b) 

and  the  constituent  relation 

Bj  =  +  Mj)  ,  (Cl.c) 

•*  -4 

where  Jj  is  an  induced  conduction  current,  Mj  is  magnetization,  and  the  relations 
Bj  =  V  x  Aj  and  Jj  =  —  a  dk^jdi  respectively  express  density  of  magnetic  flux,  Bj  ,  and 

•4  -4 

induced  current,  Jj  ,  in  terms  of  a  vector  potential,  Aj  .  We  find  from  Equations  Cl.a, 
Cl.b,  and  Cl.c  that  the  expression 


V  -  (Ai 


Hj  -  Aj 


-4 

x  Hi)  =  /io(Hj .  Mi  -  Hi- Mj)  -  <r[Aj 


Aj*  Jsi  —  Ai*  J 


sj 


(CLd) 

relates  fields  from  the  current  sources  JSi  and  JSj  . 

For  sinusoidal  current  oscillations  in  a  linear,  isotropic  medium,  the  terms  in 
Equation  Cl.d  that  come  from  magnetization  and  induced  currents  vanish  because  of 
symmetry,  so 

V  •  (Ai  *  Hj  —  Aj  x  Hi)  =  Aj  *  J  si  —  Ai*Jgj  (C2.a) 

For  confined  current  sources,  fields  vanish  at  infinity,  so  the  intergral  of  their  divergence 
over  all  space  vanishes  as  well.  Equation  C2.a  then  says  that  the  expression 


J  Aj*  J8i  d  $  =  f  Ai*  JSj  d$ 


(C2.b) 


gives  a  reciprocity  relation  between  vector  potentials  and  impressed  current  densities, 
where  intergrals  extend  over  volumes  of  respective  current  sources. 

For  filamentary  source  currents  in  circuits  i  and  j  ,  Equation  C2.b  says  that  the 
magnetic  flux  threading  circuit  i  from  unit  current  in  circuit  j  equals  the  magnetic  flux 
threading  circuit  j  from  unit  current  in  circuit  i  .  Namely,  mutual  inductance  between 
the  circuits  is  symmetric.  The  integral  expression  extends  symmetry  of  mutual  inductance 
to  volumetric  current  sources. 


Cl 


